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Abstract. A two dimensional lattice gas model with "core-softened" potential is 
investigated. Two liquid phases and density anomaly are found. The demixing 
phase transition between the two liquid phases end at a tricritical point that is also 
the terminus of a critical line. The density anomaly is shown to be related to this 
continuous line. 
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1. Introduction 

The phase behavior of single component systems as particles interacting via the so-called 
core-softened {CS) potentials is receiving a lot of attention recently. These potentials 
exhibit a repulsive core with a softening region with a shoulder or a ramp These 
models originates from the disere of constructing a simple two-body isotropic potential 
capable of describing the complicated features of systems interacting via anisotropic 
potentials. This procedure generates models that are analytically and computationally 
tractable and that one hopes are capable to retain the qualitative features of the real 
complex systems. 

The physical motivation behind these studies is the recently acknowledged 
possibility that some single component systems display coexistence between two different 
liquid phases |ni|S]-|II- This opened the discussion about the relation between presence 
of two liquid phases, the existence of thermodynamic anomalies in liquids and the form 
of the potential. The case of water is probably the most intensively studied. For 
instance, liquid water has a maximum as a function of temperature in both density and 
compressibility [Hj. It has been proposed some time ago that these anomalies might be 
associated with a critical point at the terminus of a liquid-liquid line, in the unstable 
supercooled liquid region [S], at high pressures. This hypothesis has been supported by a 
varied experimental data [H] JD] that show that thermodynamic singularities are present 
in supercooled water, around 228i^ and at atmospheric pressure. In spite of the limit of 
2357^ below which water cannot be found in the liquid phase without crystallization, two 
amorphous phases were observed at much lower temperatures JT]. There is evidence, 
although yet under test, that these two amorphous phases are related to two liquid 
phases in fluid water [12] [IE! • 

Water is not an isolated case. There are also other examples of tetrahedrally bonded 
molecular liquids such as phosphorus P [E] and amorphous silica jTHI that also are good 
candidates for having two liquid phases. Moreover, other materials such as liquid metals 
[T6] and graphite ^7) also exhibit thermodynamic anomalies. 

Acknowledging that CS potentials may engender a demixing transition between 
two liquids of different densities, a number of CS potentials were proposed to model 
the anisotropic systems described above. The first suggestion was made many years 
ago by Stell and coworkers in order to explain the isostructural solid-solid transition 
ending in a critical point fH]-Pn]- Debenedetti et al. [20] using general thermodynamic 
arguments, confirmed that a CS can lead to a coefficient of thermal expansion negative 
and consequently to density anomaly. This together with the increase of the thermal 
compressibility has been used as indications of the presence of two liquid phases [21] [22] 
which may be hidden beyond the homogeneous nucleation. The difficulty with these 
approaches is that continuous potentials usually lead to crystallization at the region 
where the critical point would be expected. 

In order to circumvent this problem, we study the effect of CS potentials in a 
lattice. Even thought the lattice is not realistic, it allows us to explore the phase space 
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in a easier way. In this work we analyze a two dimensions lattice gas with nearest- 
neighbors repulsive interactions and next-nearest-neighbors attraction. The system is 
in contact with a reservoir of temperature and particles. We show that this very simple 
system exhibits both density anomaly and two liquid phases. However, instead of having 
a critical point ending the coexistence line between the two liquid phases as one usually 
would expect, it has a tricritical point. The connection between the presence of criticality 
and the density anomaly is also shown. 

The reminder of the paper goes as follows. In sec. 2 the model is presented and 
the zero temperature phases are introduced, the mean field analysis is shown on sec. 3, 
results from simulations are discussed in sec. 4 and our findings are summarized in sec. 
5. 

2. The model and its ground state 

Our system consists of a two-dimensional square lattice with N sites. Associated to 
each site there is an occupational variable, (Tj. If the site is occupied by a molecule, 
(Tj = 1, otherwise Uj = 0. Each site interacts with its nearest neighbors with repulsive 
interactions and with its next-nearest neighbors with attractive interactions ( see Fig.l). 
Therefore the Hamiltonian of this system is given by 

h = -v,y: - ^2 E (1) 

where, < i,j > represents the sum over the nearest neighbors and << i,i » is the 
sum over the next nearest neighbors. Our system is in contact with a temperature and 
particle reservoirs. The grand potential is given by: 

^ = {n)- TS (2) 

where 7i contains the internal energy and the contribution due to the chemical potential 
/i, namely 

N 
i 

Let us now consider the ground state properties of this model. The Hamiltonian 
Eq. © allows for a number of different configurations, however, due to the lattice 
symmetry and the nature of interaction, just five of them might exhibit lower energy as 
the chemical potential is varied. They are ( see Fig. (2)) 

• Dense liquid (dl): 

<pdi = ^ = -2V,-2V2-fi. (4) 

• Uniformly diluted liquid (udl) : 

(5) 



N ^ 2 
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Figure 1. Schematic form of the interaction potential. 



• Structured diluted liquid (sdl): 

• Semi-diluted liquid (semi-dl): 

(Psemi-dl = j;^ = - Vi - V2 - -H . (7) 

• Gas (gas): 

0.a. = %^ = O. (8) 

Here (p is the grand potential per site. 

Comparing these expressions for different chemical potentials leads to the following 
zero temperature phase-diagram. For positive chemical potential, /i » \Vi\ and 
II » V2, the lower grand potential is associated with the dense liquid phase. As 
the chemical potential is reduced, the interactions between molecules become relevant. 
The first-neighbors repulsion together with the second-neighbors attraction favors the 
formation of the udl phase. At /x = — 4Vi — 2V2 there is a phase transition between the 
dense liquid phase and the uniformly diluted liquid phase. If the chemical potential is 
decreased even further, ai /j, — — 2V2 there is a transition between the uniformly diluted 
liquid phase and the gas phase. 

3. Mean-Field 

Now, let us exam the phase-diagram for nonzero temperatures employing a mean-field 
approximation. The symmetry of the different phases can be better visualized if the 
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Figure 2. (a) dense liquid phase, (&) uniformly dense liquid phase, (c) semi-diluted 
liquid phase, [d) structured diluted liquid phase, e gas and (/) fluid phase. 



Figure 3. The sub-lattices. 



square lattice is divided into four sub-lattices as illustrated in Fig. (3). In this case, the 
density of the sub-lattice a is given by: 
4 

P« = ^ E (9) 

where the sum j G a is over one of the sub-lattices a = 1,2, 3, 4. Note that the density of 
each sub- lattice varies between and 1. Using the sub-lattice division, the Hamiltonian 
given for the Eq. Q can be written as: 

^ = -EE/^'^(M)^^ (10) 
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where 



/3=1 ie/3 

is the effective chemical potential acting in an ideal sub-lattice a. Here 



= I E 



(12) 



is the chemical potential contribution due to the particles reservoir, while the second 
contribution in Eq. (|TT| is due to the interaction with other sub-lattices. The mean-field 
approximation we employ is to take the average of this last term, namely: 

4 4 



/3=i ie/3 



/9=1 



where 



P/3 



is the density of the sub-lattice (3 and where 

ea/3 = Jij, i ea, j e P 

is an interaction parameter. The mean-field Hamiltonian then becomes 



I N 



(14) 



(15) 



(16) 



/9=1 



a=l i^a \(5=1 

where the second term corrects for over-counting. It is straightforward to show that 
using Eq. (|T6|) the mean field approximation for the grand potential per site is given by 



In cosh 



(3 ^ 

13=1 



■EE 

a=l (3=1 



(17) 



- o E E ea/3P/3 + Po 
^ a=l \l3=l 

The sub-lattice density can be derived both from Eq. ()14|) and from the mean field 
grand potential by deriving it with respect to the chemical potential, namely 

I a = l,...,4 . (18) 



Po 



dpa 



In both cases, the derivation leads to; 



po 



— tanh 
2 



^ E i^c.ppf3 + Pa] 
^ /3=1 



1 

2' 



a 



1,2,3,4 



(19) 



The phase-diagram, illustrated in Fig. (4), is obtained from solving Eq. (fT^ and 
Eq. (jl7|) for different temperatures and chemical potentials. Here Vi = — 1 and V2 = 1 
are fixed. Different choices of the interaction parameters do not change the phase- 
diagram qualitatively. At high temperatures, all the sub-lattice densities are pa = 1/2 
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Figure 4. Phase Diagram for V2/IV1I = 1. The empty diamonds are continuous hnes, 
the sohd Hnes are first order transitions and the filled diamonds are the tricritical 
points. 



and the system is in the fluid phase. As the temperature is decreased at high chemical 
potential, /2 = Vl| > 2, the sub-lattices become full and the system goes from the fluid 
to the dense liquid phase. For negative chemical potentials, /2 < —2, the sub-lattices 
get empty, Pa = 0, and the system goes from the fluid to the gas phase. Between 
these two limits, 2 > /i > — 2, as the temperature is decreased two opposite sub-lattices 
become empty while the other two get full. The system goes from the fluid phase 
to the udl phase through a continuous phase transition, Tc(/i). The low temperature 
coexistence lines between the udl phase and the dl phase and between the udl and the 
gas phase are obtained by comparing the grand potentials per particle, Eq. (|T7jl . of 
these phases. The udl-dl phase boundary at /i = 2 and the udl-gas coexistence line at 
fi = —2 meet the critical line Tc{fl) at two tricritical points at (Tj = 0.924, /ij = 2) and 
(ft = 0.924, fit = -2) respectively. 



4. Monte Carlo Simulation 

The rather simple mean field approach employed in the previous session is unable to 
account for the density anomaly. For investigating the possibility of a density anomaly 
in our potential, Monte Carlo simulations in grand canonical ensemble were performed. 
The Metropolis algorithm was used to study square LxL lattice and V2/IV1I = 1. 
Different system sizes L = 10, 20, 30, 50 were investigated. Equilibration time was 
1000000 Monte Carlo time steps for each lattice site. 

The Monte Carlo simulations at fixed chemical potentials give the following results. 
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(c) id) 

Figure 5. Densities of the (a) sub-lattice 1, (6) sub-lattice 2, (c) sub-lattice 3 and {d) 
sub-lattice 4 for the lattice 20a;20 at jl = 1.2. 

At high temperatures, the density of each sub-lattice is around pa = 1/2 with sites 
randomly occupied in each sub-lattice, so the system is in the fluid phase. As the 
temperature is decreased at a fixed chemical potential, ft > 2, the density of each sub- 
lattice increases continuously to one. The system goes from the fluid phase to the dense 
liquid phase with no phase transition. If the system is cooled from the fluid phase at 
low chemical potential, p, < —2, the density decreases continuously from the fluid phase 
to the gas phase with no phase transition. As the temperature is decreased at fixed 
chemical potential, —2 < /i < 2, one finds that two opposite sub-lattices ( 1 and 3 
for example) become full while the other two sub-lattices ( 2 and 4) get empty what 
characterizes the udl phase. Fig. (5) shows that at /2 = 1.2 the density of each sub-lattice 
jumps from the fluid phase value to the udl phase value at the critical temperature 
fc{j2) = Tc/\Vi\ — 1.18 Simulations for various fixed chemical potentials, allow us to 
find the critical hne Tc{j2). 

In principle, the fluctuations observed in Fig. (5) suggest that the transition at Tc(/2) 
is continuous. This assumption is supported by the increase in the specific heat at a 
fixed chemical potential. Fig. (6) illustrates this increase at T'c(/i) = 1.18 for jj, — 1.2. 
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Figure 6. Specific iieat for tiie lattice 20x20 and ji = 1.2. 
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Figure 7. Fourtli-order Binder's cuniulant for the lattice 20a:20 and for ji = 1.2. 



The presence of a minimum in the fourth-order Binder cumulant given by 
1 - (H^) 



Vr 



(20) 



is an indication of criticahty. Fig. (7) shows Vl for /2 = 1.2 obtained from our 
simulation. The minimum at Tc(/i) = 1.1787 is a sign of the presence of criticahty. 
The hypothesis that the transition at Tc(/i) is first-order is ehminated by computing the 
energy histogram shown in Fig. (8). If the transition would be first-order, two distinct 
peaks should be present. 

In order to check the location of the coexistence lines observed by the mean- 
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Figure 8. Energy versus Monte Carlo steps, tmc for the lattice 20x20, /i = 1.2 and 
T = 1.18. The presence of only one peak characterizes a continuous transition. 




Figure 9. p.vs.ji: the first-order transitions between the gas-udl and udl-dl are 
illustrated. 



field analysis, simulations varying the chemical potential at a fixed temperature were 
performed. The results are shown in Fig. (9). For temperatures within the interval 
< T < 0.5, as the chemical potential is increased the system exhibits two first-order 
phase transitions, one between the gas and udl phase and another between the udl and 
the dense liquid phase. 

The total density of the system at fixed chemical potential is illustrated in Fig. (10). 
Differently from the mean-field result, simulations show an anomalous behavior of the 
density, instead of being monotonic with temperature, the density at positive chemical 
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Figure 10. p.vs.T: the maximum in the density for the lattice 20a;20 and fj, = 1.4 is 
shown. 




Figure 11. p.vs.T: the minimum in the density for the lattice 20a;20 and chemical 
potential jl= —1.4 is shown. 

potential increases as the temperature is decreased, it has a maximum at a temperature 
of maximum density {Ttmd{P')) and then decreases. For negative chemical potentials, 
the complementary effect is shown in Fig. (11), defining a temperature of minimum 
density {frmDifi))- 

Summarizing the results discussed above, the /i.vs.T phase-diagram is shown in 
Fig. (12). The critical line that separates the fluid from the uniformly diluted phase joint 
the phase boundaries between the uniformly diluted phase and the dense liquid and gas 
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Figure 12. Ji.vs.T phase-diagram for the lattice 20x20. The soHd hnes are the first- 
order transitions, the circles are the critical line, the filled squares are temperatures 
of maximum density (TMD), the empty squares are the temperatures of minimum 
density (TmD) and the diamonds are the tricritical points. 



phases at symmetric tricritical points at {p, = 2,T = 0.5237) and (/i = —2, T = 0.5237) 
respectively. The lines of temperature of maximum density and minimum density at 
constant chemical potential are shown. 

The p.vs.T phase diagram (see Fig. (13)) is constructed by numerically integrating 
simulations at constant pressure. It exhibits two liquid phases, a critical line and two 
tricritical points like in Fig. (12). Close to the critical line, there is a line of temperature 
of maximum density at constant pressure. 



From finite size scaling analysis |24jj, it is possible to estimate the critical 
temperature of an infinite system by the following expression: 

fo(/i)=T,(/i)(l-fa;oL-i/2) (21) 

where To(/i) is the critical temperature of the finite system, Tc(/i) is the critical 
temperature of infinite system and xq a parameter. 

Fig. (14) shows the critical temperatures of the finite systems as a function of the 
system size L. The value of Tq was obtained from two different methods: the maximum 
of the specific heat at a fixed chemical potential ft = 1.7 and the minimum of the Binder 
cumulant at the same chemical potential. The extrapolated critical temperatures are 
Tc(/2) = 0.92079 and Tc(/i) = 0.92543 respectively. The difference between these two 
values is within our error bars of our simulations. 
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Figure 14. Minimum of the Binder cumulant (dashed hne) and maximum of the 
specific heat (sohd Hnc) as a function of the system size L = 10, 20, 30, 50 for a fixed 
chemical potential ft = 1.7. 



5. Conclusions 

In this paper the phase-diagram of a two dimensional lattice gas model with competing 
interactions was investigated by mean-field analysis and Monte Carlo simulations. It 
was shown that this system exhibits two liquid phases and a line of density anomalies. 

The relation between the criticality and the density anomaly in this model goes 
as follows. The two liquid and the gas phases appear as a result of two competing 
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interactions: the softened core that favors the formation of the uniformly dilute liquid 
phase and the chemical potential that induces the system to become completely filled 
or empty. For each soft-core interaction parameter, there is a limit chemical potential 
beyond which the system is in the dense liquid phase and another negative chemical 
potential beneath which the system is in the gas phase. The double criticality arises 
from the competition between Ji and Vi. 

In systems dominated by short-range attractive forces the density increases on 
cooling. For the soft-core potential studied here similar behavior is only observed when 
the short-range repulsion becomes irrelevant ( high temperatures, /i > 2 and /i < —2). 
For 2 > /2 > —2, the soft-core repulsion prevents the density to increase to one as 
the temperature is decreased. Therefore, the same competition responsible for the 
appearance of two liquid phases leads to the density anomaly. Similar analysis can 
be made when the pressure instead of the chemical potential is kept constant. 

The presence of a critical line instead of a single critical point as one could generally 
expect jS] is not surprising. Due to the lattice structure, the udl is not one single 
phase but a region where two different phases coexist: alternating empty/full rows and 
alternating empty/full columns. These two phases become critical together with the 
dense liquid phase at the tricritical point that is also the locus the where critical line 
ends. The link between competing interactions and the presence of density anomaly 
and two liquid phases is being tested in other simple models (SEj- 
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